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Abstract 

In this paper, we model financial markets with semi-Markov volatilities and price 
covarinace and correlation swaps for this markets. Numerical evaluations of vari- 
nace, volatility, covarinace and correlations swaps with semi-Markov volatility are 
presented as well. The novelty of the paper lies in pricing of volatility swaps in 
closed form, and pricing of covariance and correlation swaps in a market with two 
risky assets. 
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Introduction 

One of the recent and new financial products are variance and voaltility swaps, which 
are useful for volatility hedging and speculation. The market for variance and volatility 
swaps has been growing, and many investment banks and other financial institutions are 
now actively quoting volatility swaps on various assets: stock indexes, currencies, as well as 
commodities. A stock's volatility is the simpliest measure of its riskiness or uncertainty. 
Formally, the volatility <jr is the annualized standard deviation of the stock's returns 
during the period of interest, where the subscript R denotes the observed or 'realized' 
volatility. Why trade volatility or variance? As mentioned in [5], 'just as stock investors 
think they know something about the direction of the stock market so we may think we 
have insight into the level of future volatility. If we think current volatility is low, for the 
right price we might want to take a position that profits if volatility increases'. 

In this paper, we model financial markets with semi-Markov volatilities and price 
covarinace and correlation swaps for this markets. Numerical evaluations of varinace, 
volatility, covarinace and correlations swapswith semi-Markov volatility are presented as 
well. 

Volatility swaps are forward contracts on future realized stock volatility, variance swaps 
are similar contract on variance, the square of the future volatility, both these instruments 
provide an easy way for investors to gain exposure to the future level of volatility. A stock's 
volatility is the simplest measure of its risk less or uncertainty. Formally, the volatility ctr 
is the annualized standard deviation of the stock's returns during the period of interest, 
where the subscript R denotes the observed or "realized" volatility. 

The easy way to trade volatility is to use volatility swaps, sometimes called realized 
volatility forward contracts, because they provide pure exposure to volatility(and only to 
volatility). 
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A stock volatility swap is a forward contract on the annualized volatility. Its payoff at 
expiration is equal to 



where o~r(S) is the realized stock volatility (quoted in annual terms) over the life of 



at is a stochastic stock volatility, K vo \ is the annualized volatility delivery price, and N is 
the notional amount of the swap in dollar per annualized volatility point. The holder of a 
volatility swap at expiration receives N dollars for every point by which the stock's realized 
volatility ctr has exceeded the volatility delivery price K vo \. The holder is swapping a fixed 
volatility K vo i for the actual (floating) future volatility <tr. We note that usually N = al, 
where a is a converting parameter such as 1 per volatility-square, and I is a long-short 
index (+1 for long and -1 for short). 

Although options market participants talk of volatility, it is variance, or volatility 
squared, that has more fundamental significance (see Demeterfi et al [llj). 

A variance swap is a forward contract on annualized variance, the square of the realized 
volatility. Its payoff at expiration is equal to 



where a^(S) is the realized stock variance (quoted in annual terms) over the life of the 
contract, 



K var is the delivery price for variance, and N is the notional amount of the swap in dollars 
per annualized volatility point squared. The holder of variance swap at expiration receives 
N dollars for every point by which the stock's realized variance a^S) has exceeded the 
variance delivery price K var . 



N{a R {S) - K vol ), 



contract 




N(a 2 R (S)-K var ), 




4 



Therefore, pricing the variance swap reduces to calculating the realized volatility 
square. 

Valuing a variance forward contract or swap is no different from valuing any other 
derivative security The value of a forward contract P on future realized variance with 
strike price K var is the expected present value of the future payoff in the risk-neutral 
world: 

P = E{e- rT (a 2 R (S)-K var )}, 

where r is the risk-free discount rate corresponding to the expiration date T, and E 
denotes the expectation. 

Thus, for calculating variance swaps we need to know only E{aj i (S)}, namely, mean 
value of the underlying variance. 

To calculate volatility swaps we need more. From Brockhaus and Long [6] approxi- 
mation (which is used the second order Taylor expansion for function ^Jx) we have (see 
also Javaheri et al [T7]): 

where V := cr^S) and 3^^3/2 is the convexity adjustment. 

Thus, to calculate volatility swaps we need both E{V} and Var{V}- 
The realised continuously sampled variance is defined in the following way: 

T 



,'idt. 



V:=Var(S):=^J a, 
Realised continuously sampled volatility is defined as follows: 



Vol(S) := ^Var{S) = W. 

Option dependent on exchange rate movements, such as those paying in a currency 
different from the underlying currency, have an exposure to movements of the correlation 



5 



between the asset and the exchange rate, this risk may be eliminated by using covariance 
swap. Variance and volatility swaps have been studied by Swishchuk |20j. The novelty 
of this paper with respect to [20] is that we calculate the volatility swap price explicitly, 
moreover we price covariance and correlation swap in a two risky assets market model. 

A covariance swap is a covariance forward contact of the underlying rates S 1 and S 2 
which payoff at expiration is equal to 



NiCovniS^S 2 )-^,. 



where K cov is a strike price, N is the notional amount, Covr(S 1 , S 2 ) is a covariance 
between two assets S 1 and S 2 . 

Logically, a correlation swap is a correlation forward contract of two underlying rates 
S 1 and S 2 which payoff at expiration is equal to: 



N{Corr R {S\S 2 )-K corrh 



where Corr(S l , S 2 ) is a realized correlation of two underlying assets S 1 and S 2 ,K corr is 
a strike price, N is the notional amount. 

Pricing covariance swap, from a theoretical point of view, is similar to pricing variance 
swaps, since 

Cov R (S\ S 2 ) = l/Ho-US'S 2 ) - o-US'/S 2 )} 

where S 1 and S 2 are given two assets, c^S) is a variance swap for underlying assets, 
Covr(S 1 , S 2 ) is a realized covariance of the two underlying assets S 1 and S 2 . 

Thus, we need to know variances for S^-S 2 and for S^-fS 2 (see Section 4.2 for details). 
Correlation Corrn^S 1 , S 2 ) is defined as follows: 
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where Cov^S 1 , S 2 ) is defined above and crf^S 11 ) in section 3.4. 

Given two assets and Sf with t e [0,T], sampled on days to = < ti < £2 < ■■■ < 

gi 

t n = T between today and maturity T, the log-return each asset is: R\ := log( /' ), % = 

S H-i 

1,2,..., ra, j = 1,2. 

Covariance and correlation can be approximated by 



Cov n (S\S 2 ) 



n 

i=l 

and 

Corr^S 1 , S 



^Var n (S^)^Var n (S 2 ) 
respectively. 

The literature devoted to the volatility derivatives is growing. We give here a short 
overview of the latest development in this area. The Non-Gaussian Ornstein-Uhlenbeck 
stochastic volatility model was used by Benth et al. [lj to study volatility and variance 
swaps. M. Broadie and A. Jain |1] evaluated price and hedging strategy for volatility 
derivatives in the Heston square root stochastic volatility model and in [5] they compare 
result from various model in order to investigate the effect of jumps and discrete sampling 
on variance and volatility swaps. Pure jump process with independent increments return 
models were used by Carr et al [7] to price derivatives written on realized variance, and 
subsequent development by Carr and Lee [8]. We also refer to Carr and Lee [9j for a good 
survey on volatility derivatives. Da Foneseca et al. [10] analyzed the influence of variance 
and covariance swap in a market by solving a portfolio optimization problem in a market 
with risky assets and volatility derivatives. Correlation swap price has been investigated 
by Bossu [2j and [3j for component of an equity index using statistical method. The paper 
|12| discusses the price of correlation risk for equity options. Pricing volatility swaps under 
Heston's model with regime-switching and pricing options under a generalized Markov- 
modulated jump-diffusion model are discussed in [13] and [14] . respectively. The paper 
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|16] considers the pricing of a range of volatility derivatives, including volatility and vari- 
ance swaps and swaptions. The pricing options on realized variance in the Heston model 
with jumps in returns and volatility is studied in [TJJj. An analytical closed-forms pric- 
ing of pseudo- variance, pseudo- volatility, pseudo-covariance and pseudo-correlation swaps 
is studied in |21j. The paper [22] investigates the behaviour and hedging of discretely 
observed volatility derivatives. 

The rest of the paper is organized as follows. Section 1 is devoted to the study of semi- 
Markov volatilites and their martingale characterization. Section 2 deals with pricing of 
variance and volatility swaps and numerical evaluation of them for markets with semi- 
Markov volatilities. Section 3 studies the pricing of covarinace and correlation swaps for 
a two risky assets in financial markets with semi-Markov volatilities. We also give here a 
numerical evaluation of these derivatives. Appendix presents a first order correction for 
a realized correlation. 

1 Martingale Representation of Semi-Markov Processes 

Let (Q, JP, (J^t)ieR+ > 5 s ) be a filtered probability space, with a right-continuous filtra- 
tion (^t)teR + and probability P. 

Let (X, 3£) be a measurable space and 

Q S M(x,B,t) := P(x,B)G x (t) for x G X,B G 3T,t G R+, (1) 

be a semi-Markov kernel. Let (x n ,T n ;n G N) be a (X x R +) SC <8> e^ + )-valued Markov 
renewal process with Qsm the associated kernel, that is 

P(x n+1 G B, r n+1 -r n <t\^ n ) = Q S Ai(x n , B, t). (2) 

Let define the process 

v t := sup{n G N : r n <t} (3) 
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that gives the number of jumps of the Markov renewal process in the time interval (0, t] 
and 

On -=T n - r n _i (4) 

that gives the sojourn time of the Markov renewal process in the n-th visited state. The 
semi-Markov process, associated with the Markov renewal process (x n , r n )„ g N, is defined 
by 

x t ■= Xu(t) for t G R + . (5) 

Associated with the semi-Markov process it is possible to define some auxiliaries processes. 
We are interested in the backward recurrence time (or life-time) process defined by 

7 (t) := t - t v ® forteR+. (6) 

The next result characterizes backward recurrence time process (cf. Swishchuk |20j). 
Proposition 1. The backward recurrence time {'j{t))t is a Markov process with generator 

Q y f(t)=f(t)+\(t)[f(0)-f(t)\, (7) 
where A(t) = G~{t) = 1 - G x {t) and Domain(Q y ) = 

Proof. Let t be the present time such that ^(t) = t, without loss of generality we can 
assume that t < n, then for T > t we have 

E i {/( 7 (T))} = E i {/( 7 (T))I ei>T } + E t {/( 7 (T))I ei < T }. (8) 

Using the properties of conditional expectation we obtain 



E*{/(7(T))} =f(T)=^ + =^—Mf(l(T))It< ei <T} 

= /(r) li + 5^)f /(r -" )G - Md "- 



(9) 
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By adding and subtracting f(t) in the integrand we get 



Mf(l(T))} = f(T)2^- + ^-r [ T (f(T - u) - f(t))G' x (u)du + f(t) Gx %.f {T \ 
then 

Mf(l(T))} - /(f) = (f(T) - f(t))=^ + ^ f\f{T ~ u) - f(t))G' x {u)du. (10) 
Now we recall the definition of the generator and using the above equation we have 

(WW = ft WMLM = m - ||[/(o) - / W ], (in 

this concludes the proof. □ 

Remark 1. ^4s well known, semi-Markov process preserve the lost- memories property only 
at transition time, then (xt)t is not Markov. However, if we consider the joint process 
(xt, 7(t))t e R +; we record at any instant the time already spent by the semi-Markov process 
in the present state, then it result that (x t , 7(t))teR + is a Markov process. 

The following result allow us to associate a martingale to any Markov process we refer 
to Swishchuk [20], Elliott and Swishchuk [15] for details. 

Proposition 2. Let (x t )tm + be a Markov process with generator Q and f G Domain (Q), 
then 

rt 

J - 



m{ := f(x t ) - f(x ) - / Qf(x s )ds (12) 

Jo 

is a zero-mean martingale with respect to & t := a{y(s); < s < t}. 

Let us evaluate the quadratic variation of this martingale (see Swishchuk [20], Elliott 
and Swishchuk |15|). 
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Proposition 3. Let (x t )tm+ be a Markov process with generator Q, f £ Domain (Q) and 
(m{)teR + its associated martingale, then 

(m f ) t := / [Qf(x s ) - 2f(x s )Qf(x s )]ds (13) 
Jo 

is the quadratic variation of . 

The next result concern the evaluation of the quadratic covariation, for details and 
proof we refer to Salvi and Swishchuk [IS]. 

Proposition 4. Let (x t )teR + be a Markov process with generator Q, f,g £ Domain (Q) 
such that fg £ Domain (Q). Denote by (wf)t 6 R + , (mf)teR + their associated martingale. 
Then 

(f(x.),g(x.)) t := f {Q{f{x s )g{x s )) - [g(x s )Qf(x s ) + f(x s )Qg(x s )]}ds (14) 

is the quadratic covariation of f and g . 

We would like to study the martingale associated to the Markov process (xt, 7(t))teR+ 
and its generator the following statement concern this task, it is a direct application of 
proposition [2j 

Lemma 1. (Swishchuk J20j) Let (xt)t&R + be a semi-Markov process with kernel Qsm 
defined in (T7]j. Then, the process 

m{ := f(x t , 7 (t)) - [ Qf(x s , 7 (s))ds (15) 
</o 

is a martingale with respect to the filtration '■= o~{x s ,u s ;0 < s < t}, where Q is the 
generator of the Markov process (x t , 7(t))teK+ given by 

Qf(x,t) = %{x,t) + 0&- [ P(x,dy)[f(y,0)-f(x,t)], (16) 
at (j x \t) Jx 

hereg x {t) = d ^. 
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The following statement follows directly from proposition Eland it allows us to evaluate 
the quadratic variation of the martingale associated with (xt,j(t))teR+- 

Lemma 2. Let (x t ) t( z^ + be a semi-Markov process with kernel Qsm, (%t, 'j(t))teR+ ^ s a 
Markov process with generator Q, f e Domain (Q) and (m{) te ^ + its associated martingale, 
then 

J\. . = f*\nf*{ 



(m f ) t := / [Qf 2 (x s , T (s)) - 2f(x s , i(s))Qf(x s , j(s))]ds (17) 
is the quadratic variation of . 

2 Variance and Volatility Swaps for Financial Markets 
with Semi- Markov Stochastic Volatilities 

Let consider a Market model with only two securities, the risk free bond and the stock. 
Let suppose that the stock price (S t )teR + satisfies the following stochastic differential 
equation 

dS t = S t (rdt + a(x t , j(t))dw t ) (18) 

where w is a standard Wiener process independent of (x, 7). We are interested in studying 
the property of the volatility a(x,j). Salvi and Swishchuk [18] have studied properties 
of volatility modulated by a Markov process here we would like to generalize their work 
to the semi-Markov case. First of all we study the second moment of the volatility, see 
Swishchuk [20J for details and proof. 

Proposition 5. Suppose that a G Domain (Q). Then 

E{a 2 {x ul {t))\^ u } = a 2 {x u , 1 {u))+ [ QE{a 2 {x s ,j{s))\^ u }ds (19) 

J u 

for any < u < t. 
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Remark 2. From proposition^ we can directly solve the equation for K{o- 2 (x t ,'j(t))\^ u } 
and we obtain 

E{a 2 (x t , 7 (t))|^J = e^a 2 (x u , 7 («)) (20) 

for any < u < t. 

In this semi-Markov modulated model we are able to evaluate high order moment of 
volatility, e.g next result concern the fourth moment. 

Proposition 6. Suppose that a 2 G Domain (Q). 

E{a 4 (x(, 7 (t))} = e*V(x,0) forteR + . (21) 

Proof. If a 2 G Domain (Q) , then from Proposition [2] we have that 

mf :=a 2 (x ul (t))-a 2 (x,0)- [ Qa 2 (x sjl (s))ds (22) 

J o 

is a zero-mean martingale with respect to & t := er{x s , 7(3); < s < t}, where Q is the 
infinitesimal generator of the Markov process (x t , 7(t))te[o,T] (cf. lemma [2]). Then cr 2 
satisfies the following stochastic differential equation 

da 2 (x t ,-f(t)) = Qa 2 (x u i(t))dt + dmf . (23) 

By applying the Ito's Lemma we obtain 

da\x ul (t)) = 2a 2 (x ul (t))da 2 (x ul (t)) + d(m° 2 ) t . (24) 

We note that, the quadrating variation of m ff2 is (see Propositon [3]) given by 

{rrf\:= [ [Qa\x s , 7 (s)) - 2a 2 (x s , 7 (s))Qa 2 (x s , -y(s))]ds. (25) 
Jo 



2.1 Pricing of Variance Swaps 
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Substituting the expression for the quadratic variation (m a ) and for a 2 , we obtain that 
cr 4 satisfies the following stochastic differential equation 

da\x ul {t)) = Qa\x ul {t))dt + 2a 2 {x ul {t))dmf . (26) 

By tacking the expectation of both side of these equation we obtain 

E{a\x t ,^{t))} = a\x,0) + [ QE{a\x s , i{s))}ds. (27) 

Jo 

Solving this differential equation we finally get 

E{a 4 (x t , 7 (t))} = e*V(x,0). (28) 

□ 

It is known that the market model with semi-Markov stochastic volatility is incom- 
plete, see Swishchuk [20]. In order to price the future contracts we will use the minimal 
martingale measure, we refer to Swishchuk [2D] for the details. Let us now focus on the 
evaluation of the price of variance and volatility Swaps. 

2.1 Pricing of Variance Swaps 

Let's start from the more straightforward variance swap. Variance swaps are forward 
contract on future realized level of variance. The payoff of a variance swap with expiration 
date T is given by 

N(a 2 R (x) - K var ) (29) 

where cr R {x) is the realized stock variance over the life of the contract 

1 f T 

a 2 R {x) := — J a 2 (x s ,-f(s))ds, (30) 



2.2 Pricing of Volatility Swaps 
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K var is the strike price for variance and N is the notional amounts of dollars per annualized 
variance point, we will assume that N — 1 just for sake of simplicity. The price of the 
variance swap is the expected present value of the payoff in the risk-neutral world 

P var {x) = E{e- rT (a 2 R (x) - K var )}. (31) 

The following result concerns the evaluation of a variance swap in this semi-Markov volatil- 
ity model. We refer to Swishchuk [20] for details and proof. 

Theorem 1. (Swishchuk 120]) The present value of a variance swap for semi-Markov 
stochastic volatility is 



yJ (e tQ a 2 (x,0)-K var )dt 



Pvar(x) = e~ rl j- J (e^a 2 (x, 0) - K mr )dt \ (32) 
where Q is the generator of (x t ,j(t))t, that is 

Q /(M ) = f ( M ) + J^- [ P(x,dy)[f(y,0)-f(x,t)]. (33) 



dV 1 ' G x (t) Jx 

2.2 Pricing of Volatility Swaps 

Volatility swaps are forward contract on future realized level of volatility. The payoff 
of a volatility swap with maturity T is given by 

N(a R (x) - K vol ) (34) 

where q-r{x) is the realized stock volatility over the life of the contract 



a R (x) := \l ^ I cr 2 (x s ,7(s))rfs, (35) 



o 



r-T 

T 

K vo i is the strike price for volatility and N is the notional amounts of dollars per annualized 
volatility point, as before we will assume that N = 1. The price of the volatility swap is 
the expected present value of the payoff in the risk-neutral world 

P vol (x) = E{e- rT (a R (x) - K vol )}. (36) 



2.2 Pricing of Volatility Swaps 
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In order to evaluate the volatility swaps we need to know the expected value of the square 
root of the variance, but unfortunately, in general we are not able to evaluate analytically 
this expected value. Then in order to obtain a close formula for the price of volatility 
swaps we have to make an approximation. Using the same approach of the Markov case 
(see also Brockhaus and Long [6 J and Javaheri at al. [17J), from the second order Taylor 
expansion we have 



Then, to evaluate the volatility swap price we have to know both expectation and variance 
of a 2 R (x). The next result gives an explicit representation of the price of a volatility swap 
approximated to the second order for this semi-Markov volatility model. 

Theorem 2. The value of a volatility swap for semi-Markov stochastic volatility is 




(37) 



Pvol (-^) 



e~ rT | yj± / T e*«cr 2 (a;,0)dt 



fo e t( ia 2 (xfi)dt) 



Var{a 2 R (x)} 




where the variance is given by 




and Q is the generator of (x t ,7(i)) t; that is 




(38) 



Proof. We have already obtained the expectation of the realized variance, 




(39) 



then it remains to prove that 



Var{a 2 R (x)} = ^ f f {e sQ [a 2 (x, 0)e tQ a 2 (x, 0)] - [e tQ a 2 (x,0)] [e sQ a 2 (x, 0)] } dsdt 



o Jo 



2.2 Pricing of Volatility Swaps 
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The variance is, from the definition, given by 



Var{a R (x)} = E{[a 2 R (x)-E{a R (x)}] 2 }, 
using the definition of realized variance, and Fubini theorem, we have 



(40) 



Var{a 2 R (x)} = E 



~jf a 2 (x tn (t))dt-^ J E{a 2 {x ul {t))}dt 



E 



2^2 



- / (a\x ul {t))-n°\xul{t))})dt 



Tj 

T r T 



(41) 



J 



E { [a 2 (x t , 7 (fj) - E{a 2 (x t , 7 (f))}] [a 2 (x s , 7 (s)) - E{a 2 (x s , T (s))}] } dsdt 



We note that the integrand is symmetric in the exchange of s and t. We can divide the 
integration on the plan in two areas above and below the graph of t=s, thanks to the 
symmetry the contribution on the two parts is the same. Then we obtain 



T r t 



Var{a R (x)} =— \ \ E{[a 2 (x t , 7 (t)) - E{a 2 (x t , 7 (i))}][a 2 (x s , j(s)) - E{a 2 (x s , i(s))}]}dsdt 



'o jo 

2 r T r* 
Wo I 



[E{a 2 (x t)7 (t))cr 2 (x s , 7 ( S ))} - E{ ( r 2 (xi, 7 (t))}E{a 2 (x s , 7 ( S ))}] dsdt. 



We would like to stress that, in this representation, the integration set is such that the 
inequality s < t holds true. Using the property of conditional expectation and proposition 
El we have 

Var{a 2 R (x)} = A £ J* [E{a 2 (x s , 7 ( S ))E{a 2 (x 4 , j{t))\& s }} - (e tQ a 2 (x,0)) (e sQ a 2 (x , 0))} dsdt. 

The process (x t)7 (t))i is Markov, then using remark [2] the conditional expected value in 
the integrand, can be expressed as 



E{a 2 (x tll (t)Ws} = e^a 2 (x s ,j(s)) =: g{x.^{ 8 )). 



(42) 



Thus, the variance becomes 



Var{a 2 R {x)} 



J>2 



T r t 



J 



[E{a 2 (x„ 7 ( S ))^(x s , 7 ( S ))} - (e tQ a 2 (x,0)) (e sQ a 2 (x,0))] dsdt. 



2.3 Numerical Evaluation of Variance and Volatility Swaps with 
Semi- Markov Volatility 
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Solving the expectation on the right hand side we obtain 

Var{*l(x)} = \[ ^ (<^%M)) - (e tQ cr 2 (x,0)) (e^ 2 (x, 0))] dsdt. 

We notice that function g evaluated in (x, 0) is simply given by 

g (x,0) =e (i - s)c V(x,0), (43) 
then substituting in the previous formula the variance finally becomes 

Var{a 2 R {x)} = ^ f f {e sQ [a 2 (x, 0)e^- s)Q a 2 (x, 0)] - [e tQ a 2 (x, 0)] [e sQ a 2 (x, 0)] } dsdt. 
T Jo Jo 

□ 

2.3 Numerical Evaluation of Variance and Volatility Swaps with 
Semi- Markov Volatility 

In the application when we attempt to evaluate the price of a variance or a volatility 
swaps we have to deal with numerical problems. The family of exponential operators 
(e*^)t involved in theorems [T] and [2] for the semi-Markov stochastic volatility model is 
usual difficult to evaluate from the numerical point of view. To solve this problem, we 
first look to the following identity 

**/(•) = E^/o, (44) 

z — ' nl 

n=0 

for any function / G Domain(Q). This identity allow us to obtain the operator {e t( ^)t at 
any order of approximation. For example, for n — 1, we obtain 

e^f{.)^{I + tQ)f{.) (45) 

where / is an identity operator. At this order of approximation we allow semi-Markov 
process to make at most one transition during the life time of contract. If we think to 
semi-Markov process as a macroeconomic factor this can be plausible. However we can 



2.3 Numerical Evaluation of Variance and Volatility Swaps with 
Semi- Markov Volatility 
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always evaluate the error in this approximation using the subsequent orders. Using the 
first order approximation the variance swap price becomes 

rT 



Pvar(x) ^e- rT |^ (I + tQ)a 2 (x,0)dt-K var ^ 
=e~ rT !^a 2 (x, 0) + ^Qa\x, 0) - K w 
Using the same approximation the volatility swap price can be expressed as 

Pvd(x) *e- rT J ± \\l + tQ)v*(x, 0)dt - — ^1 _ K 

{ V Jo 8^fi(I + tQ)<j*(x : 0)dt) 



(46) 



=e~ rT { \lo*(x, 0) + T -Qo\x, 0) - Var ^ X » _ Kvol 

Here, the variance of realized volatility is given by 

Var{a 2 R (x)} « A J jf ' {(/ + s Q)[a 2 (x, 0)(J + (t - s)Q)a 2 (x, 0)] 
- [(/ + tQ)a 2 (x, 0)] [(/ + sQ)a 2 (x, 0)] } dsdt. 



(47) 



Solving the product and keeping only the terms up to the first order in Q, we obtain 

rT ft 



Var{a 2 R (x)}^^ [ [ {sQa*(x, 0) - 2a 2 (x, 0)sQa 2 (x, 0)}dsdt 

T Jo Jo (48) 

{Qa\x, 0) - 2a 2 (x, 0)Qa 2 (x, 0)} . 



T 
3" 



Finally at the first order of approximation in Q the volatility swap price becomes 

[Qa A (x,0) -2a 2 {x,0)Qa 2 ( 
24 (a 2 (x,0) + fQa 2 (x,0)) 



V ^ 94 I rr 2 ( t. 
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3 Covariance and Correlation Swaps for a Two Risky 
Assets in Financial markets with Semi- Markov Stochas- 
tic Volatilities 

Let's consider now a market model with two risky assets and one risk free bond. Let's 
assume that the risky assets are satisfying the following stochastic differential equations 

(49) 

dSl 2) = S (2 \pf ] dt + crM(x tn (t))dwi 2) ) 
where fj,^ are deterministic functions of time, (w^) t and (w^) t are standard Wiener 
processes with quadratic covariance given by 

d[wl 1 \wl 2) ]= Pt dt, (50) 

here p t is a deterministic function and (w^) t , (wf^) t are independent of (x, 7). 

In this model it is worth to study the covariance and the correlation swaps between the 

two risky assets. 

3.1 Pricing of Covariance Swaps 

A covariance swap is a covariance forward contract on the underlying assets and 
which payoff at maturity is equal to 

N(Cov R (S^,S^)-K cov ) (51) 

where K cov is a strike reference value, iV is the notional amount and Covr(S^\ S^) is 
the realized covariance of the two assets S^ 1 ) and given by 

Cov R (S^,S^) = ipnS^lnS®] = ± £ p t a^(x tn (t))a^(x tn (t))dt. (52) 

The price of the covariance swap is the expected present value of the payoff in the risk 
neutral world 

P cov (x) = E{e- rT (Cov R (S w ,S {2) ) - K cov )}, (53) 



3.1 Pricing of Covariance Swaps 
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here we set iV = 1. The next result provides us an explicit representation of the covariance 
swap price. 

Theorem 3. The value of a covariance swap for semi-Markov stochastic volatility is 

PUx) = e~ rT jl jf Pt e tQ [*U(x, 0)a^(x, 0))dt - K cov } , (54) 



where Q is the generator of (xt,j(t))t, that is 

Qf( Xjt ) = f(x,t)+SM I P(x,dy)[f(y,0)-f(x,t)]. (55) 



dt G x (t) Jx 

Proof. To evaluate the price of covariance swap we need to know 

E{Cov R (S^\S^)} = )- [ T p t E{v il) (xt,j(t))* m (x t ,j(t))}dt. (56) 

1 Jo 

It remains to prove that 

E{a«(x*, 7 (0y (2) (^7(i))} = e tQ [a^(x,0)a^(x,0)}. (57) 

By applying the Ito's lemma we have 
d(a«(x t , 7 (t)y< 2 >(x t ,7^^ 

(58) 

+d(a^(x., 1 (-)),a^(x., 1 (-))) t . 
Using proposition H] we obtain 

d(a^\x. n (-)),a^\x. n (-))) t =Q(a^\x tn (t))a (2 \x t , 1 (t)))dt 

(59) 

-[aW(x t ,7W)Q^ (2) (^,7W) + ^ (2) (^ ) 7W)^ (1) (^,7W)]^ 

Furthermore we have 

da^(x t ^(t)) = Qa®(x u j(t))dt + dm a ^ i = 1,2. (60) 

Substituting (1591 and (1601) in equation (1581) we get 

^(Xi^W)^^^^)))^^^^^))^^^^)))^ 

(61) 

(x t , 7(t))rfm CT(2) + a (2 \x t , 7 (t))rfm ff(1) . 



3.2 Pricing of Correlation Swaps 
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Taking the expectation on both side we can rewrite the above equation as 

E{aU(x t , 1 (t)y 2 \x t , 1 (t))}=aU(x,Oy 2 \x,0) 

r* (62) 
+ / Q^ {1 \xsMs)> {2 \xsMs))}dt. 
Jo 

Solving this differential equation we obtain 

E{aM(x tn (t))<jW(x tn (t))} = e'V^MV^O)], (63) 
this conclude the proof. □ 

3.2 Pricing of Correlation Swaps 

A correlation swap is a forward contract on the correlation between the underlying 
assets S* 1 and S 2 which payoff at maturity is equal to 

N(Corr R (S\S 2 )-K corr ) (64) 

where K corr is a strike reference level, N is the notional amount and Corr^S 1 , S 2 ) is the 
realized correlation defined by 

C^,*)= , (65) 



here the realized variance is given by 

af(x) = 1 j\a^{x ul {t))fdt i = 1,2. (66) 

The price of the correlation swap is the expected present value of the payoff in the risk 
neutral world, that is 

P corr (x) = E{e- rT (Corr R (S\ S 2 ) - K corr )} (67) 

where we set iV = 1 for simplicity. Unfortunately the expected value of Corr R (S l , S 2 ) is 
not known analytically. Thus, in order to obtain an explicit formula for the correlation 
swap price, we have to make some approximation. 



3.2 Pricing of Correlation Swaps 
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3.2.1 Correlation Swap made simple 

First of all, let introduce the following notations 

X=Cov R (S\S 2 ) 

Y=af(x) (68) 
Z=af(x), 

and with the subscript we will denote the expected value of the above random vari- 
ables. Following the approach frequently used for the volatility swap, we would like to 
approximate the square root of Y and Z at the first order as follows 



Vy^Jy + y ~ y ° 



z-z 



(69) 



2v% 

The realized correlation can now be approximated by 



x 



Corr R {S\ S 2 ) « -. ^- r = -. r . (70) 

(v^+l^)(v^+f*) (l + ^)(l + ^) 

Solving the product in the denominator on the right hand side last term and keeping only 
the terms up to the first order in the increment, we have 



Corr R {S\S 2 )~ ^ oV ^° " 



I ( Y ~ Y ° + Z ~ Z ° 
2Yn 2Zn 



(71) 



l+(^ + ^) 

In what follows, we will consider only the zero order of approximation, we discuss the first 
correction on appendix. Here, we are going to approximate the realized correlation as 

Corr^S 1 ^ 2 )^^^ (72) 
v r ov 

Substituting X, Y and Z we obtain 

Corr^S 1 . S')« 1 1 / /Ha M( X „ 7 (t)y i) (x t , 7 (t))dt (73) 
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where (cf. teorem[T|), we have 



E\af (*)) = e(1 I* {a®{x ul {t))?dt\ = i T e^^^x, 0)) 2 rft, (74) 



for i = 1,2. In order to price a correlation swap we have to be able to evaluate the 
expectation of both side of equation (173|) . the expectation of the right hand side becomes 

E {a^(x t , 7 (t)V (2) (x t , 7 (t))} = e tQ a^(x, 0)a^(x, 0). (75) 
We can summarize the previous result in the following statement. 
Theorem 4. The value of a correlation swap for semi-Markov stochastic volatility is 
_ rT { / Q T p t e^[a ( - 1 \x,0)a^(x, 0)]dt \ 

Pcorr\%) ~ 6 \ , / Kcorr / i C^) 

I V lo ^(aW^O)) 2 ^ ^ e t( 3(a( 2 )(x,0)) 2 rft J 
where Q is the generator of {x t ^{t))t, that is 

Qf( x ,t) = f (x,t) + 0$- [ P(x,dy)[f(y,0)-f(x,t)]. (77) 
at (j x \t) Jx 

4 Numerical Evaluation of Covariance and Correlation 
Swaps with Semi-Markov Stochastic Volatility 

In order to obtain a more handy expression for the price of covariance and correlation 
swaps to use in the application, we will introduce here an approximation for the family of 
operator (e t ®)t- Following the approach used for the variance and volatility case, we are 
going to approximate the operators at the first order in Q as 

e^f{.)^{I + tQ)f{-). (78) 

Using this approximation the covariance swap price becomes 

Pcov{x) ~e- rT [f£pt( J + tQ)[a^{x, 0)a^{x, 0)]dt - K ( 



r /•" r T i (79) 

|a (1) (x,0)a (2) (a;,0)y p t dt + Q[a^(x,0)a ( - 2) (x,0)) J tp t dt - K co A . 
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The same approximation allow us to express the correlation swap price as 

p t (I + tQ)[aW(x,0)aW(x,0)}dt 



■0*0 



-rT 



-K r . 



=e 



-rT 



yfjjjXl + tQ)(<rM(x, OWdtyJffil + tQ)(*W(x, 0))*dt 

gWjx, 0)q( 2 \x, 0) f[ptdt + Q[aW(x, 0)a^(x, 0)] £ tp t dt 
y/(<rW(x, 0))2 + f Q(aW(x, 0))^(a^(x, 0))* + ZQ{aW(x, 0))< 



(80) 
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Appendix 



A Correlation Swaps: First Order Correction 

We would like to obtain an approximation for the realized correlation between two 
risky assets bea 

Corr R (S\ S 2 ) = , C °" ii(S '; ^ (81) 



In section 13.2.11 we have already obtained the following approximated expression 

x 

— .X 

1 - 



Corr R (S\S 2 )- ~ X 



Y-Y , Z-Z 



21q 2Zq 



52) 



where 



F=4 1)2 (x) (83) 
Z=af(x), 

and with the pedix we have denoted the expected values. We have already evaluated 
the expectation of the zero order approximation, now we would like to evaluate the first 
order. 

Substituting X, Y and Z in equation (|82|) we obtain 

Corr R (S\S 2 ) « 1 1 / p iff < 1 '(i 1 , 7 («))ff (2, (i„7(*))* 



1 



T p it 7W(x t ,7(t))a (2) (a:t,7W)* 



2T2(E{4 1)2 (x)})3/2(E{ag )2 (x)})3/2 Jo " ( 8 4) 
x|E{4 2)2 (x)}^ T [(aW(x s , 7 ( S ))) 2 -E{(a( 1 )(x s , 7 ( S ))) 2 }R' 

+ E{4 1)2 (x)} /" T [(a( 2 )(a; u , 7 ( M ))) 2 -E{(a( 2 )(a; u , 7 (n))) 2 }]dnL 



26 



where 

EWo*(*)} =E|i jT(aW(x t , 7 (*))) 2 *} = - jTe*(a«(s,0)) 2 ^ (85) 

for i = 1,2. We have to evaluate the expectation of the right hand side of equation 
We already calculated the expectation of the first term, which is the zero order 
approximation for the realized correlation. Then we will focus now on the other terms. 
First of all, let rewrite them as follows 



p t *M(x tn (t))*W(x t ,j(t)) (E{4 2) (x)}[(a^(x sn (s))f-M{(a^(x sn (s))f}) 
o Jo v (86) 

+ E{4 1)2 (x)}[(a( 2 )(x s ,7( S ))) 2 -E{(a( 2 )(x s ,7( S ))) 2 }])M, 



we have four different contributions in the integrals, the expectation of the terms 

f / T p t( r( 1 )(xi,7(t))a (2) (x i ,7W)E{ ( T« 2 (x s ,7(s))}E{a^ 2 (x s ,7(s))}M (87) 
Jo Jo 

for i — 1, 2, can be evaluate using theorem [3j Then, in order to evaluate the expectation 
of the approximated realized correlation, it only remains to calculate 

E {So So P '° r(1)(x '' lit))a(2) ^)) aW2 l{s))dsdA t = l,2 (88) 

To this end, let's first divide the range of integration in two intervals as follows 



E Uo Jo P^&'^yKxtnWaV { Xail {s))d8& 

+ f / T p i( T( 1 )(x i ,7(t))^ (2) (^ ) 7W)^ (t)2 (a; s ,7(s))M 
Jo Jt 



39) 



for i = 1, 2. We notice that the first integral set is such that t > s while the second has 
t < s. We can now use the property of conditional expectation to obtain 



+ 



E<| / / pM{^ 1 \x t Mt))° {2 \xtMt)WsW i) {x Sll {s))dsdt 

(90) 

p,a( 1 )(x 4 , 7 (t))a^(x t ,7(t))E{a« 2 (x s ,7(s))|^}^l. 
o Jt ) 



Jo 

T pT 
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We notice that (xt,j(t))t is a Markov process then using the Markov property, we can 
express the conditional expectations as 

E{^ 1 \x t , 1 {t))a^{x hl {t))\^ s ) = e^W^x.nWaWix,,^)) =■ %„7(<0) 
for t > s, and 

E{a^(x s ,j( S ))\^ t } = e^a^\x ul {t)) =: g®{x ul {t)) 

for s > t. Therefore, the first term of eq. (l90l) can be expressed as 

e(/ / M(:r s , 7 (s))a (i)2 (x s ,7(s))dsdA = / / p<e sQ [/i(x, 0)a {i) \x, 0)}dsdt, (91) 
I Jo Jo J Jo Jo 



while the second as 

fT r T 



E 



o Jt 

T pT 



Pt o( l \x t Mt)W 2 \xtMt))g^(xtMWsdt\ 

p t e tQ [<x (1) (x, 0)a (2) (a;, 0)g (i) {x, 0)]dsdt. 



(92) 



'0 Jt 

Now, we can evaluate the functions h and g at x obtaining 

h(x, 0) = e {t ~ s)Q (x, 0)<7< 2 > (x, 0)] (93) 

and 

2«(x,0) = e (s -' )Q [ ( x« 2 (x,0)]. (94) 

We can summarize the previous result in the following statement which gives the correla- 
tion swap price up to the first order of approximation. 

Theorem 5. The value of the correlation swap for a semi- Markov volatility is 

Pcorr(x) = e~ rT (E{Corr R (S\ S 2 )} - K corr ) (95) 
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where the realized correlation can be approximated by 

2/ T p t e tQ a (1 \x,0)o-W(x,0)dt 



E{Corr R (S\S 2 )} 



^// T e«3((r( 1 )(a;, 0)) 2 dty/ e«3(aW(x, 0)) 2 dt 



Jo Ptilo eSQ i etQ [ (t(1) ( x > ) a{2) ( x > °)] a(1)2 ( x > °) } ds + It e tQ W (1) fa 0)aW {x, 0)e"« [<7« 2 (x, 0)] }du)dt 

3/2 / ~ TT72 (96) 



2 (fe'Q^M^.O)) 2 ^ (/ T e*«( ( 7( 2 )(x,0)) 2 ^ 
Jo T A(/o e sQ {e«? (x, 0)a^ (x, 0)}a^ 2 (x, 0)}ds + jf e< Q {a^ (x, 0)a^ (x, 0)e u ® [a( 2 ) 2 (x, 0)] }du)dt 



2 (/^^((tW^.O)) 2 ^) 172 (/ T e«3(<T( 2 )(x,0)) 2 dt) 3/2 
here Q is the generator of the Markov process (x t ,^(t))t given by 

Qffa t) = £(x,t) + ^&J x Pfa dy) [f(y, 0) - f(x, t)} . (97) 
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